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Abstract : Swimming, i.e., being able to advance in the absence of external forces 
by performing cyclic shape changes, is particularly demanding at low Reynolds num- 
bers which is the regime of interest for micro-organisms and micro-robots. We focus 
on self-propelled stokesian robots composed of assemblies of balls and we prove that 
the presence of a wall has an effect on their motility. To rest on what has been done 
in [I] for such systems swimming on R 3 , we demonstrate that a controllable swimmer 
remains controllable in a half space whereas the reachable set of a non fully controllable 
one is affected by the presence of a wall. 
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1 Introduction 

Swimming at low Reynolds number is now a well established topic of research which 
probably dates back to the pioneering work of Taylor (25] who explains how a micro- 
organism can swim without inertia. Later on, Pureed [19 formalized the so-called 
"scallop theorem" which states that, due to the reversibility of the viscous flow, a 
reciprocal deformation of the body cannot lead to a displacement of the swimmer. 
However, this obstruction can be circumvented using many swimming strategies |19| . 
Swimmers can be distinguished with respect to their ability to change their shape or to 
impose rotational motions of some parts of their body in order to create viscous friction 
forces on the fluid, and produce by reaction, the propulsion. 

Many applications are concerned by this problem as for example, the conception of 
medical micro devices. The book by J. P. Sauvage [21] presents a lot of engine models 
adapted for tiny devices while the design and fabrication of such engines have been 
recently investigated by e.g. B. Watson, J. Friend, and L. Yeo [26 . As an example, 
let us quote the toroidal swimmer, first introduced by Purcell [19] and which has been 
subsequently improved by A.M Leshansky and O. Kenneth [12], Y. Or and M. Murray 
|18| . A. Najafi and R. Zargar [17 a among others. 

The strategy for swimming consists in a cyclic deformation of body with a non- 
reciprocal motion. The first swimmer prototype belonging to this class is the three link 
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swimmer also designed by Purcell |19] . More recently, R. Golestanian and A. Ajdari 
|lL)j introduced the Three-sphere swimmer which is geometically simpler and allows for 
exact calculations of motion and speed [2] , or even explicit in some asymptotic regimes 

In the continuation of [2], F. Alouges, A. DeSimone, L. Heltai, A. Lefebvre, and B. 
Merlet [T] showed that the trajectory of the Three-sphere swimmer is governed by a 
differential equation whose control functions correspond to the rate of changing shape. 
The swimming capability of the device now is recast in terms of a control problem to 
which classical results apply. 

Of particular importance for applications is the issue of the influence of any bound- 
ary on the effective swimming capabilities of micro-devices or real micro-swimmers. 
Indeed, boundaries clearly affect the hydrodynamics and may have an influence on the 
swimmer's capabilities. In that direction, an biological study of Rothschild [20] claimed 
for instance that spermatozoids tend to accumulate on walls. More recently, H. Winet, 
G. S. Bernstein, and J. Head [27] proved this related boundary effect for the sperm of 
humans which evolves in a narrow channel. Swimming in a geometrically confined envi- 
ronment then became a subject of major interest, in particular to model this attraction 
phenomenon (see [23], MM)- D - J Smith, E. A. Gaffney and J. R. Blake [24 have 
described the motion of a stylized bacterium propelled by a single flagellum and they 
show that the attraction by the wall is effective. Later, H. Shum, E.A Gaffney, and J. 
Smith [22] investigated to which extent this attraction effect is impacted by a change 
in swimmer's morphology. 

On a more theoretical side, other approaches provide results that show an attraction 
effect by the wall. A. P. Berke and P. Allison [3J, modelling the swimmer with a simple 
dipole, put in evidence an attraction due to the presence of the wall. Y. Or and M. 
Murray [IS] derived the swimmer dynamics near a wall for three various swimmers, but 
with unvarying shapes. The case of a changing shape swimmer has been studied by 
R. Zargar and A. Najafi |2§j , where the dynamics of the Three-sphere swimmer in the 
presence of a wall is given. However, some fundamental symmetry are not satisfied in 
their swimmer's motion equation. 

The aim of this paper is to attack the same problem (the influence of a plane wall in 
the motion of the swimmer) by means of control theory. Several recent works present a 
controllability results for a self-propelled micro-swimmers in a space without boundary, 
as example, let us quote the paper of J. Loheac, J. F. Scheid and M. Tucsnak p3] and the 
study of J. Loheac and A. Munnier [13] made of the spherical swimmer in the whole space 
(see also |6J for the same kind of results in a perfect fluid). Furthermore, F. Alouges, A. 
DeSimone, L. Heltai, A. Lefebvre, and B. Merlet pQ deal with the controllability on R 3 
for the Three-sphere swimmer and others specific swimmers. The question that we want 
to address now is whether the presence of the plane wall modifies the controllability 
results. We here prove two results in that direction. Namely, considering the fully 
controllable Four sphere swimmer proposed in pQ , we show that in the half space, the 
swimmer remains fully controllable, while a Three-sphere swimmer enriches its reachable 
set, at least generically which seems at first sight contradictory with earlier results. 
Indeed, although previous works show an attraction from the boundary, the set of 
reachable points could be of higher dimension. In other words, if the dynamics is 
somehow more constrained due to the presence of the wall, the set of points that the 
swimmer may reach could be larger than what it was without the wall. 
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The rest of this paper is organized as follow. In Section [2] we describe the two 
model swimmers to which our analytical and numerical tools are later applied. Section [3] 
presents the main controllability results associated with the introduced swimmers. In 
Section [4] we show that swimming is indeed an affine control problem without drift by 
using a similar approach than F. Alouges, A. DeSimone, L. Heltai, A. Lefebvre, and 
B. Merlet in pQ. The controllability result is proved in Section [5] for the Four-sphere 
swimmer and in Section [6] for the Three-sphere swimmer. Concluding remarks are given 
in Section 

2 Swimmers 

We carry on the study of specific swimmers that were considered in [I] in R 3 . In order 
to fix notation, the wall is modeled by the plane W = {(x, y, z) G M 3 s. t. y = 0}, and 
the swimmers, which consist of N spheres (-Bj)i=i..iv of radii a connected by thin jacks, 
are assumed to move in the half space R3. = {(x, y, z) G R 3 s. t. y > 0}. As in [1 , the 
viscous resistance associated with the jacks is neglected and the fluid is thus assumed 
to fill the whole set R 3 ^ \ uf =l Bi. The state of the swimmer is described by two sets of 
variables : 

• the shape variables, denoted by £ (here in R^ -1 or R^), which define the lengths 
of the jacks. A stroke consists in changing the lengths of these jacks in a periodic 
manner ; 

• the position variables, denoted by p G R 3 x 5*0(3), which define swimmer's 
position and orientation in the half-space. 

In what follows, we call S C R for a suitable M G N the set of admissible states 
(£, p) that we assume to be a connected nonempty smooth submanifold of R A/ . We 
thereafter focus on two swimmers that have been considered in the literature, the Three- 
sphere swimmer (see [16] . [2J, pQ) and the Four sphere swimmer (see [I]). It turns out 
that this latter is easier to understand that the former, and we therefore start with it. 

2.1 The Four-sphere swimmer 

We consider a regular tetrahedron (Si, S2, S3, S4) with center O G R+. The swimmer 

consists on four balls linked by four arms of fixed directions OSj which are able to 
elongate and shrink (in a referential associated to the swimmer). The four ball cluster 
is completely described by the list of parameters (£,p) = (£1, . . . , £4, c, TV) G S = 

(y^oo) 4 x R.^ x SO(3) C (^JJa, oo) A x R 6 (a rotation TZ G 50(3) is uniquely 
characterized by its 3 dimensional rotaion vector). It is known (see pQ) that the Four 
sphere swimmer is controllable in R 3 . This means that it is able to move to any point 
and with any orientation under the constraint of being self-propelled, and when the 
surrounding flow is dominated by the viscosity. This swimmer is depicted in Fig. [TJ 
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Figure 1: The Four-sphere swimmer. 
2.2 The Three-sphere swimmer 

This swimmer is composed of three aligned spheres as shown in Fig. [2] We assume that 
at t = the swimmer starts in the vertical half-plane H = {(x, y, 0) S R 3 s. t. z = 
0,y > 0}, it is clear from the symmetry of the problem that the swimmer stays in H 
for all time, for whatever deformation of its arms it may carry out. We characterize 
swimmer's position and orientation in H by the coordinates (c, 9) G R 2 x [0, 2ir], where 
c G H is the position of one of the three spheres, and 9 is the angle between the swimmer 
and the x— axis. Therefore, in that case, the swimmer is completely described by the 
vector (£, p) = (&, c, 9) G S = (2a, oo) 2 x H x [0, 2vr) C (2a, oo) 2 x R 2 x R/2vrZ. In 
the three dimensional space R 3 (when there is no boundary), it is obvious by symmetry 
that the angle 9 cannot change in time, and thus this swimmer is not fully controllable. 
One of the main contributions of this paper is to understand the modifications of this 
behavior due to the presence of the plane wall. 




Figure 2: The Three-sphere swimmer. 
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3 The main results 

Consider any of the swimmers described in the previous Sections, and assume it is self- 
propelled in a three dimensional half space viscous flow modeled by Stokes equations. 
In this paper, we will establish that both swimmers are locally fully controllable almost 
everywhere in S. By this we mean the precise following statements. 



Theorem 3.1 Consider the Four-sphere swimmer described in Section 2.2, and assume 
it is self-propelled in a three dimensional viscous flow modeled by Stokes equations in 
the half space R+. Then for almost any initial configuration (£\p 4 ) G S, any final 
configuration (£^,p^) in a suitable neighborhood of (£ l ,p*) and any final time T > 0, 
there exists a stroke £ G W 1,oo ([0, T]), satisfying £(0) = £ l and £(T) = and such that 
if the self-propelled swimmer starts in position p l with the shape at time t = 0, it 
ends at position p? and shape at time t = T by changing its shape along £(t). 



Theorem 3.2 Consider the Three-sphere swimmer described in Section 2.1, and as- 
sume it is self-propelled in a three dimensional viscous flow modeled by Stokes equations 
in the half space Rj_. Then for almost any initial configuration (£\p l ) G S such that 
p % G H, any final configuration (£^,p^) in a suitable neighborhood of (£*,p*) with 
pf € H and any final time T > 0, there exists a stroke £ G W 1 ' oo ([0, T]), satisfying 
£(0) = £ l and £(T) = and such that if the self-propelled swimmer starts in position 
p* with the shape at time t = 0, it ends at position pf and shape at time t = T by 
changing its shape along £(t) and staying in H for all time t G [0, T], 



Remark 3.1 The sense of "almost every initial configuration" can be further precised 
as everywhere outside a (possibly empty) analytic manifold of codimension 1. 

The proof of the controllability of the Four-sphere swimmer is given in Section [5] 
whereas Section [6] is devoted to demonstrate devoted to demonstrate Theorem 13.21 



4 Mathematical setting of the problem 

As for their 3D counterparts, the equation of motion of both swimmers take the form 
of an affine control problem without drift. In this section, we detail the derivation of 
this system. 



4.1 Modelization of the fluid 

The flow takes place at low Reynolds number and we assume that inertia of both the 
swimmer and the fluid is negligible. As a consequence, denoting by £1 = \jf =l Bi the 
space occupied by the swimmer, the flow in Rjj_ \ satisfies the (static) Stokes equation 

' -fiAu + Vp = in R^O, 

div u = in R^ \ Q,, 
< —cm = f on d£l, (1) 
u = on 9R^, 
u — > at oo. 
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Here, we have denoted by <x = /i(Vu + V*u) — pld the Cauchy stress tensor, n is 
the unit normal to dfl pointing outward to the swimmer. We also set 

V = {u G V'(K 3 + \ n, R 3 ) | Vu G L 2 (R 3 \ ft), = G L 2 (R 3 \ ft)} . 

V 1 + l r l 

It is well known that V is a Hilbert space when endowed with the norm (and the 
associated scalar product) 

Hull 2 ; := I |Vu| 2 . 



We also assume that f G H 1//2 (<9ft) in order to obtain a unique solution (u,p) to 
the problem QmVx L 2 (M 3 f \ ft) which can be expressed in terms of the associated 
Green's function (obtained by the method of "images", see (5]) as 



u(r) = / K(r,s)f(s)ds, (2) 

Jan 

where the matricial Green function K = j=i 2,3 is given by 

K(r, r ) = G(r - r ) + Ki(r, r ) + K 2 (r, r ) + K 3 (r, r ) , (3) 

the four functions G, Ki, K2 and K3 being respectively the Stokeslet 

1 /Id rigir" 
87r/x \ |r| |r| 



G(r) = ^(£ + ^) (4) 



and the three "images" 

1 /Id r' <g> r' 
8-Kfi \ \r'\ |r 



Kl (r,r ) = -— ( — + ^3-) , (5) 



^i( r ' r o) = ^o 2 (l-Mfe-^) , (6) 



*3,y(r, r ) = -j^-yo (1 - 2S j2 ) ( ^U, - jlfa + ^S j2 - 3 ^ ) . (7) 



Here ro = (xo,yo,zo) and r' = r — fo, where fo = (xq, —yo,zo) stands for the "image" 
of ro, that is to say, the point symmetric to ro with respect to the wall. 

Let B be the sphere of radius 1 centered at the origin. We identify the boundary 
of the domain occupied by the swimmer, <9ft, with (dB) N and we represent by fj G 
H~ 1 /' 2 (dB) the distribution offeree on the sphere Bi. Correspondingly, Uj G H l l 2 (dB) 
stands for the velocity distribution on the sphere Bi (and of the fluid due to non-slip 
contact). 

Following we denote by 7(£ >p ) the Neumann-to-Dirichlet map 

%P)= ^" 1/2 ^ (8) 

(fi,...,fjv) h-» (ui,...,u N ) 

where we have denoted by H.^ 1 / 2 the space 

(iJ±l/2(aB))iV It ig well known that the 
map 7(£ iP ) is a one to one mapping onto while its inverse is continuous. 
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Using ([2]), we can express u« (i = 1,2,3) by 

N 

Vr G dB, Uj(r) = Y/ K(xj + ar, x,- + as)fj(s)ds 



TV 

^(fj,K(xj + ar,Xj + a-))as 



(9) 



where ( v )flB stands for the duality (H" 1 / 2 (dB), H 1 / 2 (dB)) . 

Proposition 4.1 The mapping (£, p) 1-4 77g iP ) is analytic from S into C(l~L~ l l 2 ,1-L l l 2 ) . 
Furthermore, 77^ jP ) is an isomorphism for every (£, p) G S, and the mapping (£, p) 1-4 
is a/so analytic. 

Proof: The proof is identical to the one given in [fj, replacing the the Stokeslet by the 
Green kernel K which is also analytic outside its singularity. □ 



Remark 4.1 As the direct consequence, the mapping 77| !P ) and its inverse depends 
analytically on a. 

4.2 Equation of motion 

In this section, we use the self-propulsion assumption in order to express the dynamics 
of the swimmer as an affine control system without drift. 

Proposition 4.2 There exists a family of vectorfields Fj G T5, such that the state of 
the swimmer is described by the following ODE, 

1( p)5<- do) 

Proof: 

This equation of motion is by now classical in this context (see [I], [2], [8] or |15j). 
Let us recall the principle of its derivation. 

At any time t, the swimmer occupies a domain (we therefore denote by Oo the 
domain occupied by the swimmer at time t = 0). We also define the map <!> which 
associates to the points of 80,q x S, the current point in <9Qt, 

* : dn xs -4 on t 

(x ,£,p) 1-4 x t . 

When inertia is negligible, self-propulsion of the swimmer implies that the total viscous 
force and torque exerted by the surrounding fluid on the swimmer vanish i.e., 

(12) 
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From the linearity of the Neumann-to-Dirichlet map, we deduce that the system 
(12) reads as a linear system which depends on p and £. By inverting it, we get p 
linearly in terms of Assuming that £ £ R fc for some k £ N, we thus obtain 

k 

P = £ W *&P)&> (13) 

i=i 



canonical basis. 



which becomes (10) when we call Fj := ( ], where Ej is the i-th vector of the 



Wi 



□ 



Let us recall some notations which are used to study the controllability of such 
systems of ODE (see for instance |11|). 

Let F and G be two vector fields defined on a smooth finite dimensional manifold 
Ai. The Lie bracket of F and G is the vector field defined at any point X £ Ai 
by [F, G]{X) := (F • V) G(X) - (G • V) F(X). For a family of vector fields I on M, 
Lie(J-) denotes the Lie algebra generated by T . Namely, this is the smallest algebra 
- defined by the Lie bracket operation - which contains T (therefore T C Lie(J-) and 
for any two vectorfields F £ Lie(F) and G £ Lie(F), the Lie bracket [F, G] £ Lie(J-)). 
Eventually, for any point X £ Ai, Liex(F) denotes the set of all tangent vectors V{X) 
with V in Lie{F). It follows that Liex{J~) is a linear subspace of TxAi and is hence 
finite-dimensional. 

Lie brackets and Lie algebras play a prominent role in finite dimensional control 
theory. Indeed, we recall Chow's theorem: 

Theorem 4.3 (Chow |7J) Let Ai be a connected nonempty manifold. Let us assume 
that J- = (Fi) 7 ^ =1 , a family of vector fields on Ai, is such that Fi £ C°°(Ai, TAA) , Vi £ 
{!,••• ,m}. 

Let us also assume that 

Liex(F) = T X (M) , VX £ M . 



Then, for every (X° , X 1 ) £ M*M, andforeveryT > 0, thereexistsu £ L°°([0, T); R m ) 
such that the solution of the Cauchy problem, 

m 

x=j: Ui f(x), (14) 

X(0) =x°, 
is defined on [0,T] and satisfies X(T) = X . 



The theorem 4.3 is a global controllability result, we also recall the one which gives 
a small-time local controllability. 

Theorem 4.4 ([7J, p. 135) Let 0, be an nonempty open subset of R n , that T = 
{Fi)™ =1 , a family of vector fields, such that Fi £ C°°(J7, R n ) , Vz £ {1, • • • , m} . 
Let X e such that 

Hex AT) =R n - 
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Then, for every e > 0, there exists a real number rj > such that, for every (A^X 1 ) £ 
{X s. t. \\X — X e \\ < rj} 2 , there exists a bounded measurable function u : [0,e] — > R n 
such that the solution of the Cauchy problem 

m 

* = 5>*w. (15) 

X(0) = X° , 

is defined on [0, e] and satisfies X(e) = X . 

When the vector fields are furthermore analytic (and the manifold M is also analytic) 
one also has the Hermann-Nagano Theorem of which we will make an important use in 
the theoretical study of the controllability for our model swimmers. 

Theorem 4.5 (Hermann- Nagano [11]) Let A4 be an analytic manifold, and T a 
family of analytic vector fields on Ai. Then 

1. each orbit of T is an analytic submanifold of M, and 

2. if N is an orbit of J 7 , then the tangent space of M at X is given by Liex(J~). In 
particular, the dimension of Liex^J 7 ) is constant as X varies over N. 

In our context, the family of vector fields is given by T = (Fj)i<j<fc which are defined 
on the manifold M = S, and the controls Ui are given by the rate of shape changes £j. 
In view of the preceding theorems, the controllability question of our model swimmers 



raised by Theorems 3.1 and 3.2 relies on the dimension of the Lie algebra generated by 
the vectorfields (Fj)i<j<^ which define the dynamics of the swimmer. In particular they 
are direct consequences of the following Lemma. 



Lemma 4.6 For almost every point (in the sense of remark 3.1) (£,p) € S, the Lie 



algebra generated by the vectorfields ^i)\<%<h a t (£>p) equal to T^ p )5. 

The proof of this lemma is developed until the rest of the paper. Several tools are 
used in order to characterize this dimension among which we mainly use asymptotic 
behavior and symbolic computations. As we shall see, although the theory is clear, the 
explicit computation (or at least asymptotic expressions) is by no means obvious and 
requires a lot of care. In particular, before using symbolic calculations, a rigorous proof 
of the expansion, together with a careful control of the remainders in the expressions 
allowed us to go further. 



5 The Four-sphere swimmer 



In this section, we give the proof of the controllability result stated in Theorem 3.1 



Proof: The argument of the proof is based on the fact that K given by ^ satisfies 

K(r,r') = G(r-r') + oQ , (16) 

where r = (x, y, z) and r' = (x 1 ', y' , z') are two points of R^_, and G is the Green function 
of the Stokes problem in the whole space R 3 , namely the Stokeslet, defined by (B. 



9 



As a consequence, we obtain that the Neumann to Dirichlet map given by Q satisfies 
for a swimmer of shape £ at position p = (p x ,p y ,p z ,lZ) E R 3 x SO3 



%p) =T? + 



(17) 



where is the Neumann-to-Dirichlet map associated to the Green function G . 
The system (12) now reads 



7? 



-1 



+ O 



Xf x 



7? 




((d p <D)p + dx t = , 

(9 p $)p + (d^)t) dxt =0. 



(18) 



Consequently, the ODE (10) becomes 



(19) 



Py 



(20) 



where (F^)j=i 5 ... ^ are the vector fields obtained in the case of the whole space R 3 . 
In other words, we obtain the convergence 

F(£,p)=F°(0 + O 

and also for all its derivatives to any order. 

It has been proved in [T] that dim Lie^(F°) = 10 at all admissible shape showing 
the global controllability in the whole space of the underlying swimmer. We thus obtain 
that for p y sufficiently large 

dim Lie (€iP) (F) = 10, (21) 



and therefore due to the analyticity of the vector fields (Fj)j=i 5 ... (21) holds in a 
dense subset of S. This shows that the system satisfies the full rank condition almost 
everywhere in S and proves Lemma 4.6 in this context, and thus Theorem |3.1| by a 
simple application of Chow's theorem. □ 



The preceding proof can be generalized to any swimmer for which the Lie algebra 
satisfies the full rank condition in R 3 . We now turn to an example for which this is not 
the case, namely the Three-sphere swimmer of Najafi Golestanian (TB]. Indeed, when 
there is no boundary, this swimmer is constrained to move along its axis of symmetry. 
The purpose of the next section is to understand to which extent this is still the case 
when there is a flat boundary. 



6 The Three-sphere swimmer 



This section details the proof of Theorem 3.2 It is organized in several subsections 



each of them focusing on a particular step of the proof. In the subsection |6.1| by 
introducing some notations, we recall the expression of the equation of motion of the 
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Three-sphere swimmer. Subsection 6.2 deals with the special symmetry which have to 
be verify by the vector fields of the motion equation. From this symmetry properties, 
we deduce the reachable set of the particular case where the swimmer is perpendicular 
to the wall. In the subsection |6.3| we give an expansion of the Neumann- To-Dirichlet 
mapping associated to the Three-sphere swimmer and its inverse, in the case where 
the radius of the sphere o is small enough and the distance of the arm is sufficiently 



large. In subsection 6.4 we deduce from this previous approximation an expansion of 



the motion equation of the swimmer, for a sufficiently small. Finally, the subsection 6.5 
presents some formal calculations of the vectors fields of the motion equation and their 
Lie brackets which leads to obtain, almost everywhere, the dimension of its Lie algebra. 

6.1 Equation of motion for the Three-sphere swimmer 



From Section |2.2[ we know that the swimmer's position is parameterized by the vector 
(x, y, 9) where (x, y) is the coordinate of the center of -B 2 as depicted in Fig. [2] and 9 is 
the angle between the swimmer and the x— axis. We recall that £ := (£1,^2) stands for 
the lengths of both arms of the swimmer. 



The motion equation (10) thus reads 



d_ 
dt 



( 6 \ 

X 

V 

\0 J 



Fi(£,x,y,0Ki + F 2 (£,x,y,0)6. 



(22) 



Notice that, from translational invariance of the problem, both Fx and F2 actually do 
not depend on x. 

In what follows, we denote by 



d 



((,v,o) 



dim Lie 



(Fi,F 2 



the dimension of the Lie algebra Lieuyfi\(F\) F2) C M 5 at (£, y, 1 
and F2 are independent one to another and never vanish, that 

2 < d (ty,o) < 5 • 

6.2 Symmetry properties of the vector fields 

Proposition 6.1 Let S be the 5x5 matrix defined by 

( 1 \ 

1 

0-10 

1 

V -1 / 

Then one has for all £ = (£1, £21 x, y,9) £ S 

F 1 (£i ) £ 2 ,y,0) = SF 2 (£ 2 ,£ 1 ,y,27r-0) 
and similarly for the Lie bracket 

[Fx , F 2 ] (£ 1; £ 2 , y, 9) = S [F 2 , Fi] (£ 2 , £1, y, 2tt - 9) . 



It is clear, since Fi 
(23) 



(24) 
(25) 
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Proof: Although the plane breaks the 3D axisymmetry along the swimmer's axis, we 
can still make use of the symmetry with respect to the vertical plane that passes through 
the center of the first sphere -B 2 . A swimmer with position (x,y,9) and shape (£i,<!; 2 ) 
is transformed to one at position (x,y,2ir — 0) and shape (^d) ( see Fig. Making 
use of the fact that corresponding solutions to Stokes equations are symmetric one to 
another, we easily get the proposition. 









f Xi \ 














X^V 1-6 






v 1 J V x ~ 3 J 







wall 



Figure 3: The plane symmetry which links the situation at (£1, £ 2 > y, ff) with those at 
(£2; £i) x, y, 2ir — 9). In both cases, solutions to Stokes flow are also symmetric one to 
another. 

Eventually, one deduces the Lie bracket symmetries by applying the former symme- 
tries on the vectorfields themselves. An easy recurrence shows that the same identities 
hold for any Lie bracket of any order of the vectorfields Fi and F2. In particular one 
has for instance 



[Fi , [Fi , F 2 ] ] (£1 , & , y ,9) = S [F 2 , [P 2 , Fi]] (6, £1 , V, 2tt 



(26) 



□ 



As a direct consequence of proposition |6.1[ we deduce that the fourth coordinate of 
the Lie bracket [Fi,F 2 ] vanishes at (£,£, y, 0) and at (£, £, y,ir). 

Proposition 6.2 Let T be the 5x5 matrix defined by 

( 1 \ 

10 

0-10 

1 

V -1 J 
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Then one has for all £ = (£1, £ 2 , x, y, 9) £ S, and i = 1, 2 



Fi(£i,£ 2 ,y,0) =TFi(£i,6,y,7r-0) 



(27) 



and similarly for the Lie bracket 



[Fi,F 2 ](£i,&,y,0) =T[Fi,F 2 ](^,e 2 ,!/,7r-e). 



(28) 



Proof: The two identities readily come from the symmetry which transforms a swimmer 
with position (x, y, 0) and a shape (£1, £2) to one at position (x, y, 0) with the same shape 
(see Fig. 



Figure 4: The plane symmetry which links the situation at (£1, £2, a?, y, 9) with those at 
(£l,£ 2 ,a;,y,7r — 6*). In both cases, solutions to Stokes flow are also symmetric one to 
another. 

Eventually, one deduces the Lie bracket symmetries by applying the former symme- 
tries on the vectorfields themselves. An easy recurrence shows that the same identities 
hold for any Lie bracket of any order of the vectorfields Fi and F 2 . □ 

As a result, in the case where 9 = ±5, we get the dimension of the Lie algebra of 
the vector field Fi and F 2 . 

Corollary 6.3 The dimension of the Lie algebra Lie^ ^ !2/>7r / 2 )(Fi, F 2 ) is less than or 



Proof: We deduce from the preceding proposition that for i = 1,2 and j = 3,5, 
F^(£i, £2, y, ±7r/2) = 0. This simply means that a swimmer starting in the vertical 
position cannot change its angle 9 and its abscissa x by changing the size of its arms. 




wall 



equal to 3. 
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As a matter of fact, the same holds true for any Lie bracket of Fi and F2 at any order, 
and we can deduce from this that 

since any vector of the Lie algebra Lie^ £ 2 ^ „■ /2) (F^ , F2) has a vanishing third and fifth 
component. 

Moreover, by using the argument introduced in Section [5] we get, for almost every 
y, that the dimension of the Lie algebra is almost equal to the one without boundary 
(see [I]). □ 



Remark 6.1 We denote, for all x 6 R, the set 

:= {(£l,&,z ,y,±|) s. t. (6,6) G (2a,oo) 2 y > oj 

which corresponds to the case where the swimmer is perpendicular to the wall and A x 
the set of states where the the dimension of the Lie algebra generated by Fi and F2 is 
equal to two, i.e., 

A x := j((£i,&,a;,y,±^)s. t.d^ 2jy ^/ 2 ) = 2j . 
5y using the property of analyticity \4-l\ Ax is a finite union, 

A * = U \{Ci,^,x,y,±^)s. G (2a, oo) 2 ) . 

We deduce, that for all x, j\f x \A x defines a set of three-dimensional orbits strictly in- 
cluded in the manifold S. 



Furthermore, the proof of the corollary 6.3 can be applied to all generic positions 



then, it implies that the dimension of the Lie algebra is almost equal to 3, almost 
everywhere, i.e., 

6.3 Approximation for small spheres and large distances 

For the general case {9 7^ 7r/2), the preceding computation is not sufficient to conclude. 
In order to proceed, we make an expansion of the vectorfields and their Lie brackets 
with respect to a (the radius of the balls) near 0. 

This part is devoted to the proof of the expansion of the Neumann to Dirichlet map 



(|34j) together with its inverse (35) at large arms' lengths. Let us first define for all 
G {1,2,3} 2 , the linear map %j as 

Ti,j : H-V 2 {dB) -> H l / 2 {dB) 

fj 1 — y / K(xj + a-, Xj + as) fj(s) ds . 

JdB 
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We recall that the Green kernel K writes (following (|3|) as 

K(r, r') = G(r - r') + Ki(r, r') + K 2 (r, r') + K 3 (r, r') 



where G is the Stokeslet (see Q) and each kernel is given by the corresponding coun- 
terpart in ([3]). Eventually, we call 7~° the Neumann to Dirichlet map associated to 
G 

f / G(a(- -s))f(s)ds. 

JdB 

Proposition 6.4 Let € {1,2,3} 2 . We have the following expansions, valid for 

a <C 1: 

• if i j then 

Tij = K(x ij x i )(f i) Id)aB + Ri (29) 
where ||Ri||£(h- 1 /2 j h 1 /2) = O (a) , 

• otherwise 

3 

Ti,i = T° + ^K fc (x i ,x i )(f i ,Id) aB + R 2 (30) 
fc=i 

w/iere H^H^^-i^^i^) = O (a) . 
Proof: Let G {1,2,3} 2 be such that i / j, and fj G H- 1 l 2 {dB). We define 

VredB, u;(r) := (77,f,)(r) = / K(x; + ar, Xi + as)f,(s)ds , (31) 

./as 

and 

Vj(r) = Uj(r) - K(xj,x,) / f,(s)cfa= / (K(xj + ar,x,- + as) - K(xj,x,-)) f,(s)ds. 

JdB JOB 

Our aim is to estimate the H l / 2 {dB) norm of Vj. Butj^] 

H v i||.ffV2(aB) - ll v i||fri(B)> 

and since K(x,y) is a smooth function in the neighborhood of x = x$ and y = Xj, one 
has Vr, s & B 

|K(xj + ar,xj + as) -K(xj,Xj)| = O(a) , (32) 
and for the gradients in both r and s 

|V r K(xj + ar, Xj + as)\ = O (a) , 
| V s K(xj + ar, Xj + as) | = (a) , 
| V r V s K(xi + ar, Xj + as) | = O (a 2 ) . 



Here and in the sequel, we use the definition for the H 1/2 {dB) norm: 

II v IIhi/2 (ss) = min ||w|| H i (s) 

v ' wsh" 1 (b,r 3 ),w=v on as 
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Therefore, we obtain Vr £ B 

|v,(r)| < ||K(xj + ar,Xj + a- 



K ( x M x i)ll^|ll f ill H -^ 



< 0(^11^11^1, 



and similarly 



|V P Vi(r)| < ||V r (K(x i + ar J x i + aO)|| if i||f i || Jf _i 



< 0(a)\\ti 



This enables us to estimate the H 2 norm of v« on dB 

llv,:ll__i._. < llv, ; " 



L3(fl) 



+ ||Vvi 



< 0{a 



which proves (29) 



In order to prove (30), we use the decomposition ([3| where none of the kernels 
(Kj)j = i 2,3 is singular. Therefore Vr G dB 

Uj(r) := (7i,ifi)(r) = / K(x; + ar, x, + as)fj(s)ds 

= / G(a(r - s))fi(s)ds + / (Ki + K 2 + K 3 )( Xi + ar, x, + as)f;(s)ds 

= T\+ / (K 1 + K 2 + K 3 )(x 4 + ar,x J + as)f J (s)ds. 

•/SB 

We finish as before, having remarked that for / = 1,2,3 

Ki(xi + ar, Xj + as) = K/(xj, Xj) + O (a) . (33) 

□ 



Proposition 6.5 For every f £ % - 1 / 2 , 

3 

(r x f) 4 (r) = T\ + ^K J (x i ,x i )(f i ,Id) SB + ^K(x J ,x i )(f J) Id) aB + 7^(f), (34) 



1=1 



with 



ith \\Ki 11^(^-1/2^1/2) = 0(a). 
Proof: 

For all i & 1,2, 3, and all r £ dB 



(T x i)Ar) := / K(xj + ar,Xj + as) fj(s)ds + ^ / K(xj + ar, x,- + as) f}(s)ds 
JdB i^j JdB 



and the result follows from the application of (29) and (30) of Proposition 6.4 



□ 
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Proposition 6.6 In the regime a <C 1, one /ias for every u £ T-L 1 ^ 2 , 

{T- l u).= (r°)- 1 (u l -ELiK fc (x l ,x i )((r°)- 1 u J ,id) aB 



Xj , X ,-)((T°)- 1 u J -,Id) SB + 7^(u) 



(35) 



with \\lZi\\ £^1/2^-1/2 ^ = O (a 3 
Proof: We recall that 



T°:i?~^(&B) -»■ m{8B) 

f ^ / G(a(--s))f(s)ds, 

and define for / = 1,2,3 the operators 

Si:H~^{dB) -> H^{dB) 

/ Ki(xj,Xj)f(s)ds, 



f i ^ 



and eventually 



Sij : H~2(dB) -»• Hi(dB) 



f ' ^ 



/ K(xj,Xj)f(s)c?s. 



That these operators are continuous operators from H~?(dB) into Hz(dB) is classical. 
We hereafter are only interested into the estimation of their norms, and more precisely 
the way they depend on a, <5 and y in the limit a — ?► 0. Notice that since the kernel G 
is homogeneous of degree -1, one has 



HT°|| 



£(tt-Va,tti/2) 



o 



and 



As far as 5; is concerned, we get that (since |K/(xj,Xj)| = O (1)) 



0(a) 



and similarly 



Oil) , 



0(1) 



(36) 

(37) 
(38) 



When a — > this enables us to invert (34) leading to (p5b. 



□ 



6.4 Self-propulsion 

We now use the fact that the spheres are non-deformable and may only move following a 
rigid body motion. In other words, the velocity of each point r of the i— sphere expresses 
as a sum of a translation and a rotation as 

Ui(r) = u Ti +u Ri (r) , (39) 
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where is constant on dB while Ui^(r) = u; x ar for a suitable angular velocity u)i 
(remember that all quantities are expressed on the unit sphere dB). This is of peculiar 
importance for the computation of the total force and the total torque, which, due to 
self-propulsion, should vanish which implies 



,)B 



dB 



0. 



(40) 



Plugging (39) in (40) and using (35) leads to 



]T / (TV 1 fu Ti +U Ri -^K fc (x i ,X i )((r°)- 1 (uT l + U a ),Id)eB N ) - 

^E K ( x - x i)((^°) _1 (uT J +u Rj ),Id)dB ] j = O (a 3 ) ||u| 



( ( ^ (41) 

It is well known that both translations and rotations are eigenfunctions of the Dirich- 
let to Neumann map of the three dimensional Stokes operator outside a sphere. Namely 



A T u Ti and IT ) u Ri = X R u R 



It is well-known that At = 3 -^ t leading in particular to the celebrated Stokes formula 



-1 



while \ R = 3fia. We also remark that due to J dB u^ds = , we have 



JdB V 



-1 



u^ds = 0. 



We therefore obtain 



6TTfia^2 [ u Ti ~ 6717/0 ^ K fc (xi,Xi)u Ti - 67r^a^ K(x i ,x J )u Ti = O (a 3 ) ||u| 



k=l 



(42) 

We now compute the torque with respect to the center xi of the first ball B\ . Self- 
propulsion of the swimmer implies that this torque vanishes: 

0= / arxfi(r)+/ (x2- Xl +ar) xf 2 (r)+ f (x 3 -xi+ar) xf 3 (r) = Ii+I 2 +I 3 , (43) 
Job Job JdB 

( cos 9 

where the quantities Ii, I 2 and I 3 are respectively given below. Calling eg = sin# 

V 
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the direction of the swimmer 
Ii 



= / (xi-x 2 + or) xfi(r) = / (free + or) x (T x ) 1 u 1 
JdB JdB 

= / (-^ee + ar) x (T^) -1 u T1 -6vr/iaVKi(xi,xi)u Tl 

JdB \ fri 

— 67r//ay^K(xi,Xj)uy 7 - + O (a 2 ) ||u|| J 
= — Qirixa^ie.Q x I u^ 1 — 67TjLta^K;(xi,xi)uy 1 — 6717x0 y^K(xi,x ? -)ur,- j + 0^a 3 ) ||u| 

Similarly, we get, 

I 2 = a / r x f 2 (r) = o / r x (Tx) -1 ^ 

= a/ r x (7 -0 ) -1 u T2 + - 67r/io VKi(x 2 ,x 2 )u T2 

V ^ 

—6717^0 K(x 2 , x^u^,- + O (a 2 ) ||u| 
= a / rx(T )- 1 (u i?1 ) + o(a 4 )||u||. 
But since (T°) 1 uri = A^uri = Ar^i x ar, we have 

a r x (TV 1 ( U R2) = a 2 X R / rx(wixr) dr 
JdB JdB 



8vr 



-jia uj\ 



This leads to 



Correspondingly, 

I3 = / (x 3 - x 2 + or) x f 3 (r) 



8tt 



fia 3 LOi + O ^a 4 ) ||u| 



67r/ia^ 2 e x u T3 - 6717x0^] K z (x 3 , x 3 )u T3 - 6717m K(x 3 , Xj)utj + 0(a 3 ) ||u| 

1=1 m 



Denoting by A the matrix 



' An A i2 A i3 N 

A 2 i A 22 A 23 
y A 3 i A 32 A 33 J 



(44) 
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where for i = 1,2,3 

3 

Ajj = Id - 6irfia K z ( x i > x 
i=i 

and for i,j = 1,2,3 with i ^ j 

Ay = — 67r^oK(xj, Xj) 



(45) 



(46) 



and S the matrix 



Id Id Id 

-£ie e x +£ 2 e0X 



we can rewrite the self propulsion assumption ( 42 ) , ( 43 ) as (notice that angular velocities 
being involved of higher order disappear) 



SA 



U T2 
\ U T 3 



O a 2 ) llu 



(47) 



We end up by expressing uy^UT^u^ and u\, 0)2,(^3 in terms of x,y, 6, £1 and £2- 
But, since is the velocity of the center of the ball Bi, one has 





' X 


- 6 cos(#) + Hi sin(fl) \ 




/ 




UTi = 


V 


- |i sin(0) - Hi cos(#) 


, U T 2 = 




S) 




V 


/ 




V 





and 



Similarly 



UT 3 



( x + 6 cos(0) - 56 sin(0) 
y + 6 sin(6») + H2 cos(0) 




/ 




We rewrite these formulas as 



with 



' UTj N 




UT 2 


l-'(i) 


V U ^3 , 






'id -a4 


T = 


Id 







(48) 



where 



- sin ( 
cos 9 I and 



V 



U 



/ -ee 

V e e 
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Plugging (48) into (47) leads to the motion equation 



x 



(SA + R) 



V + U£ = 



(49) 



V 



V 



where the residual matrices have a norm which is estimated as 

IIRII = 0(a 2 ^ 



6.5 Dimension of Lie algebra under the small spheres hypothesis 



Rewriting from (44 ), (45 ) and (46 ) A = Id+aAi, we can expand in power series of a the 



solution of (49). This enables us to write an expansion (still in a) of the two vectorfields 
Fi and F2. To this end, we have used the software MAPLE to symbolically compute 
those expressions and the Lie brackets [Fi,F 2 ], [Fi, [Fi, F 2 ]], and [Fi, [Fi, F 2 ]]. Writing 
the vectorfields in components as 



Fi(£i,&,V, 



1 





Ff + O (a 2 ) 
Ff + 0{a 2 ) 
V Ff + 0(a 2 ) J 



\ ( 

F 2 (6,6,2/,0) := 




1 



F| + O (a 2 ) 
n + O (a 2 ) 
V Fl + 0(a 2 ) j 



(50) 



we find, after having furthermore expanded the abovementioned components in power 
series of -, 

y 



1 



a 



3 a 



cos(#) + -cos(0)Kt{£ 1 ,b,0) + (sin(0) cos(#) (6 + 26)) 



3 w 6 
a 



16y 2 
a 



384y 3 
1 



(cos(0)iY 2 3 (6,6,0)) + — 4 (smWcoBfflKittu&e^+O 



a 



3 a 



cos(#) - - cos(0)l^(6,6, -0) + tttt (sm(0) cos(0) (26 + 6)) 



3— V-y g— v-y-x^-^x, . 1%2 



384y J 



(cos(0)£T 2 3 (6,6,-#)) + FTq 4 (sm(0)cos(0)#f(6,6,-0)) +0 



Ff = -sin(0) + -8^(0)1^(6,6, 

3 O 



512y 
3 a 



W sin(0)K|(6,6^) - ^(6,6,^ + 
isin(0) - ^(0)1^(6,6,-0) _ _^_^4 (6 6 ^ 



g^sin^) 008(0)1^(6,6,0) 

3 a 
64y 2 



u " co S (e)Ki(Z l ,&,6) + / " 



512y 



9a 



sin(0) 008(0)^(6, £1, -^) + ^4 cos(0)^ 2 & (6,6, -0) + O 



a 
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In those expressions the remaining functions are respectively given by 



j (6 2 + 66 + 6 2 ) 66 (6 + 6) 

iff (£, 0) = -210fi cos(6<) 2 + 12 cos(#) 4 £ 2 + 184£ 2 + 24 cos(#) 2 66 

-3266 - 6cos(0) 4 66 - 92^ 2 + 105£ 2 cos(fl) 2 - 6cos(0) 4 6 2 
Kl(tO) = m - c \ - m (12 cos(^) 4 g 2 5 + 24^ 5 cos(fl) 4 - 168g 2 5 cos(fl) 2 - 336^ 5 cos(g) 2 

+112£ 2 5 + 72&£ - 176£ 2 $ ~ 136^e 2 2 + 224^ 5 - 156^2 cos(#) 2 
-2A^l cos(fl) 2 - 240£ 4 6 cos(fl) 2 + 48^6 - 1566^2 cos(#) 2 
-24cos(#) 4 6 2 6 3 + 9cos(#) 4 M! - 2ie i 3 cos(^)% 2 ) 
K 4 ( t (&i 2 -ffii-d + 2gfe + 2#) 

j (6 2 + 66 + 6 2 ) 66 (6 + 6) 

K 4 (£,0) = 6cos(^) 2 6+3cos(^) 2 6-46-2e 2 

if 4 (6#) = -132£ 2 cos(0) 2 + 6cos(0) 4 £ 2 + 56£? + 12cos(0) 2 66 

-1666 - 3cos(fl) 4 66 - 28£f + 66£fcos(#) 2 - 3cos(0) 4 £f 

*1(6*) = (512C 2 + 51 4fe + 512CD ( " 210 C ° S(g)4$2 " 42Q ^ C ° S(g)4 + 232$l C ° S(g)2 
+24cos(0) 6 6 5 - 64£ 2 5 + 12cos(0) 6 e 2 5 - 9666" - 646 2 £ 2 3 
-128^1 + 104^^1 cos(6>) 2 - 56^1 cos(6) 2 
-96^ 2 + 2166 ?2 cos(0) 2 - 66cos(0) 4 e 2 Cf 
-318£ 4 cos(fl) 4 6 - 2406 3 cos(0) 4 ^ 2 - 24 cos(fl) 6 6 2 6 3 
-21 cos(#) 6 6 3 £ 2 + 464^ 5 cos(#) 2 - 128^^1 
+264£ 4 £ 2 cos(#) 2 - 204 cos(#) 4 66 4 + 9 cos(#) 6 6£ 2 4 ) , 



= -86 -4^ + 2 cos(6>) 2 6 + cos(#) 2 £ 2 



772 T c c X ^ ( 20 cos (^) 2 6 4 - 40 cos(fl) 2 £ 4 - 4 cos(#) 4 6 4 + 8^ cos(0) 4 
(4 2 + 66 + K\) v 

-406 3 6 - 86 2 6 2 + 32C 3 6 - 166 4 + 32£ 2 3 cos(^) 2 6 - 7cos(#) 4 6 3 6 

+ cos(^) 4 6 2 6 2 - 40 cos(#) 2 6 3 6 + 8C 3 cos(^) 4 6 + 32£ 4 - 8£&? cos(#) 2 ) . 

As one can see, the use of a software for symbolic computation seems unavoidable. 
Subsequently, we get the expansion of the Lie bracket [Fi, F 2 ](6, 6> V, ®) 



( 



[F l ,F 2 ](ty,9) :-- 








\ 



[F 1 ,F 2 ] 3 + 0(a 2 ) 
[F!,F 2 ] 4 + 0(a 2 ) 
V [F 1 ,F 2 ] 5 + 0(a 2 ) J 



(51) 
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where the components are given by the following expressions 
r F pi a PfW m (g + 2 ^ 2 + &l + + 6 4 ) 

[Fi ' F2]s - "3 cosW wrwm 

27 a acos(.) sin(^ 2 " 4 + 8 ) ^ + ^ 



512 y 4 w w fe 2 + 66 + 6 2 

+0 



[Fi,F 2 ] 4 = --sin(0) 



3 w (6 + 6) 2 e 2 2 6 2 

27a 2 66 (cos(g) 4 - 4 cos(fl) 2 + 8) (-g 2 + g) 

+ 5i2^ cosW (6 2 + 66 + 6 2 ) 

+ o ( 4 



[Fi,F 2 ] 5 = 



81a gi6(6+6)(co S W 4 -4cos(g) 2 + 8) 

4 COS^J y^iy - h C 



512 y 4 



(6 2 + 66 + e 2 ) 



Notice that since the two first coordinates of Fi and F 2 are constant, the corresponding 
first coordinates of the Lie bracket vanish. Similarly, the asymptotic expansion for the 
second order Lie bracket [Fi, [Fi, F 2 ]] (£, y, 6) reads 



/ 



[Fi, [Fi, F 2 ]] (£, y, 








\ 



[F 1 ,[F 1 ,F 2 ]] 3 + 0(a 2 ) 
[F 1 ,[F 1 ,F 2 ]] 4 + 0(a 2 ) 
V [F 1 ,[F 1 ,F 2 ]] 5 + 0(a 2 ) J 



where 



[Fi,[Fi,F 2 ]] 3 



[Fi,[F 1; F 2 ]], 



3 cosW c?(a+6) 3 



27a 
512 y 4 



cos(#) sin(0)L 3 (£,0) + 



2a . ^ 6 (3^ + 3^6 + CI) 
smlfcM — - — ^ ^ — — 

3 e 3 (a+6) 3 



cos(#) 2 L 3 (6#) + 



[F 1 ,[F 1 ,F 2 ]] 5 = 



27a 
512 

^cos(0)L 4 (*;,0) + O 



There, L 3 and L4 are respectively given by 

(8 - 4cos(#) 2 + cos(fl) 4 ) (2£ 2 - 66 - £!) 



U{i,9) = 



(6 2 + 66 + 6 2 ) 2 

6 3 (8-4cos(g) 2 + cos(g) 4 ) (26+6) 

(6 2 + 66 + 6 2 ) 2 



(52) 
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Eventually, the expansion of the vector field [F2, [Fi, F2]] is given by 



where 



[F 2 ,[Fi,F 2 ]] (Si,&,y,0) :-- 








[F 2 ,[F 1 ,F 2 ]] 3 + 0(a 2 ) 
[F 2 ,[F 1 ,F 2 ]] 4 + 0(a 2 ) 
V [F 2 ,[F 1 ,F 2 ]] 5 + 0(a 2 ) J 



[F 2 ,[Fi,F 2 ]] 3 = 



[F 2 ,[Fi,F 2 ]] 4 = 



— a cost 8) — - — 5 

3 (el(a+6) 3 ) 



27a 
512y z 



acos(0)sin(0)L 3 (£2,£i,-0) + 



2a sin(0) 6 (^ + 3 ^ + 3 ^) 



cos 



[F 2 ,[Fi,F 2 ]] E 



3 fe 3 (6 + 6) 3 

27a 
~512y 4 

-^4Cos(0)L 4 (6,a,-^) + O 



(0)%(&,a,-0) + o( 



a 

Vy 5 

a 

y 5 



(53) 



We now can compute an expansion of det (Fi,F 2 , [Fi,F 2 ], [Fi, [Fi,F 2 ]] , [F 2 , [Fi,F 2 ]]) 
which if it does not vanish implies the local controllability of our model swimmer. It 
can be readily checked that we have 



det (Fi,F 2) [Fi,F 2 ], [Fx, [Fi,F 2 ]] , [F 2 , [Fi,F 2 ]]) = 

[Fi.Fala [F l5 [F 1? F 2 ]] 3 [F 2 , [F l5 F 2 ]] 3 

[F 1; F 2 ] 4 [Fi,[Fi,F 2 ]] 4 [F 2 ,[F!,F 2 ]] 4 

[F 1; F 2 ] 5 [F 1 ,[F 1 ,F 2 ]] 5 [F 2 ,[F 1; F 2 ]] 5 



81a 3 (6-6) 
131072 y 9 



sin6>(cos0) 2 #(t#) + O 



1 

.Tlo 



(54) 



with, 



(6^ 6 + 27^6 + 50^^2 2 + 55gej + 50^i 2 g 2 4 + 27fr£f + 6£ 2 6 ) 

(6 + 6)fe 2 + £i6 + £ 2 ) 2 66 

x (64 - 64 cos(#) 2 + 32 cos(#) 4 - 8 cos(#) 6 + cos(6<) 8 ) 



It is easily seen that R never vanishes. Therefore, the previous determinant has 
a non- vanishing first coefficient (in —A which does not vanish for £1 / £2 and 8 4 

{0, §,vr, ^}. Since it is an analytic function of (£,i,£,2,y,9) we deduce that it does not 
vanish except at most on a negligible set. This is sufficient to conclude that 



almost everywhere and the local controllability of the Three-sphere swimmer around 
such points. 
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Remark 6.2 Quite strikingly, when £1 = £ 2 the first term of the expansion vanishes 
and one has to go one step further. We find in that case 

det (Fi, F 2 , [F l5 F 2 ], [F 1; [F u F 2 ]] , [F 2 , [F X ,F 2 ]]) (£, f , 2/, 0) = Tfo #)4o + ° (4l 



where 



T(£,y,< 



945 
"524288 



a 3 sin(#) 2 cos(#) 2 £ (cos(#) 4 + 8-4 cos(#) : 



This coefficient does not vanish unless 8 ^ {0, f ,vr, 



The case # = or n. We already know from symmetry that when 8 = \ or 8 = 
~, one has d(£ 1( £ 2 y0) < 3. Therefore, it remains to understand the case 8 = (or 
7r by symmetry). The preceding computation does not allow us to conclude about 
the dimension of the Lie algebra at such points. Indeed, the 2 first coefficients of 
the expansion of the determinant vanish, and it might well be the case at all orders. 
Nevertheless, in that case, we can expand the subdeterminant 



A 



Fi,[Fi,F 2 ]] 3 
[F 1} F 2 ] 5 [FijFi.FaHs 



[Fi,F 2 ] 3 



in order to obtain informations. Indeed, one gets 



A 



45 
512 y 



;a 2 (£l-6Kl#'(0 + O 



with, 



(6+6) 2 £ 2 2 fe 2 + £ 2 + £i6) 



,^1 + ueid + 16^ 2 3 + i9e 3 c 2 2 + 12^6 + 3^ 



As the direct consequence, we get that the dimension of the Lie algebra, c% )2/l o) > 4, 
for almost every (£,?/) G (R + ) 3 . 

□ 



This finishes the proof of Lemma 4.6 and thus of Theorem 3.2 



Remark 6.3 As usual, it is possible to pass from local to global controllability on each 



of the connected components where the determinant given by (54) does not vanish. More 
precisely, let A := j(£>p) s. t.d^ p ^ < 4z\, we define by <S(£ )P ) the connected component 
of the subset S\A which contains (£, p). Applying Chow's Theorem 4-3 on S((,p), gives 
that for every initial configuration (£\p*), any final configuration (£t , pf ) in <S(£ lP ), 
and any final time T > 0, there exists a stroke £ £ W 1,co ([0, T]), satisfying £(0) = £* 
and £(T) = ^ and such that the self-propelled swimmer starting in position p l with the 
shape £* at time t = 0, ends at position pf and shape ^ at time t = T by changing its 
shape along and staying in <S(£ )P ) for all time t G [0,T]. In other words, S(£ p) is 
exactly equal to the orbit of the point (£,p). 
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7 Conclusion 



The aim of the present paper was to examine how the controllability of low Reynolds 
number artificial swimmers is affected by the presence of a plane boundary on the fluid. 
The systems are those classically studied in the literature (see [T] for instance) but are 
usually not confined. This is the first in-depth control study of how the presence of the 
plane wall affects the reachable set of a peculiar micro-swimmer. 

Firstly, the Theorem 3.1 shows that the controllability on the whole space implies 
the controllability in the half space. Although the proof is applied on the Four-sphere 
swimmer, it is based on general arguments which can be appropriate for any finite 
dimensional linear control systems. 

Secondly, the Theorem 3.2 deals with the controllability of the Three-sphere swim- 
mer in the presence of the plane wall. We prove that, at least for this example, the 
hydrodynamics perturbation due to the wall surprisingly makes the swimmer more con- 
trollable. This result is not in contradiction with the several scientific studies which 
show that the wall seems to attract the swimmer (see [21]]) [21]) [22]) [2]) [I])- Although, 
the Theorem 13.21 leads to the fact that the wall contributes to increase the swimmer's 
reachable set, we can conjecture that some of them are easier to reach than others. 

The quantitative approach to this question together than the complete understand- 
ing of the situation in view of controllability of the underlying systems is far beyond 
reach and thus still under progress as is, in another direction, the consideration of more 
complex situations like, e.g. rough or non planar wall. This is the purpose of ongoing 
work. 
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